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The process that prevents the deposition of cooled gas in cooling flows must rely on 
feedback in order to maintain gas with short cooling times, while preventing the bulk 
of the gas from cooling to low temperatures. The primary candidate for the feedback 
mechanism is the accretion of cooled and cooling gas by an active galactic nucleus 
(AGN). Despite some difficulties with this model, the high incidence of central radio 
sources in cooling flows and the common occurrence of radio lobe cavities, together, 
support this view. The Bondi accretion rate for the intracluster gas onto the AGN 
depends on the gas properties only through its specific entropy and that is governed 
directly by competition between heating and cooling. This provides a viable link for 
the feedback process. It is argued that the mass accreted between outbursts by the 
central AGN is only sensitive to the mass of the black hole and the gas temperature. 
Bondi accretion by an AGN leads to a simple expression for outburst energy that can 
be tested against observations. 



1. Introduction — the Significance of AGN 
Heating 

A variety of mechanisms have been proposed to explain 
why large quantities of gas do not cool to low temper- 
atures at the centers of cooling flows (many discussed 
elsewhere in these proceedings). While there are many 
candidate heating mechanisms, in general, short central 
cooli ng times, as little as 10^ - 10^ y in ma ny clusters 
(e.g. lDavid et al.ll2001t pTakizawa et alJl2003|) and an or- 
der of magnitude shor ter still in elliptical galaxies (e.g. 
iMathews fc Brighentil |2003). are very difficult to main- 
tain without a process involving feedback. In the absence 
of feedback, any heating process will either overwhelm, 
or be overwhelmed by radiative cooling in a few cool- 
ing times. Adding to this the observation that outbursts 
from central AGN's have an impact on the the hot gas 
through the creation of cavities around radio lobes (e.g. 
i Blanton et"al]l200ltlFabian et al.l200(]l:lMcNamara et all 
|2QQ3) in clusters, groups and isolated elliptical galaxies, 
and the prevalence of radio sources at the centers of cool- 
ing flows (jB urns. ■1990'!. suggests that AGN feedback is 
the prime candidate as the process that prevents mass 
deposition in cooling flows. 

This is not to dismiss t he difficulties faced by the 
AGN f eedback model Ce.g. iBrighenti fc Mathew^ l2003t 
iFabian et al...2001t and elsewhere in these proceedings). 
However, it has been well demonstrated that a mod- 
est amount of thermal conduction could make a signifi- 
cant contribution to the energy budget of the regions in 
clusters where the cooling time is shorter than the age. 
Thermal conduction could, therefore, alleviate the ma- 
jor problem of th e AGN feedback model in clu sters, but 
as also shown by ' Zakamska fc NaravanI l)2003f) , thermal 
conduction cannot prevent the deposition of cooled gas 
in all cluster cooling flows. The steep temperature de- 
pendence of the thermal conductivity (k ~ T^/^) means 



that thermal conduction is less effective in groups, and 
less still in isolated elliptical galaxies. Furthermore, the 
region in clusters where the cooling time is shorter than 
the Hubble time is relatively small, so that the remainder 
of the cluster can act as a thermal reservoir to make up 
for radiative losses from the center. This is not the case 
in isolated elliptical galaxiess, where radiative losses are 
significant throughout the system. Thus, while thermal 
conduction may be important for the energy budget of 
cluster cooling flows, it is unlikely to have much impact 
in cooler systems. 

For the clusters with radio cavities studied to date, sin- 
gle AGN outbursts have very little impact on the cluster 
as a whole. However, the opposite is true for isolated 
elliptical galaxies. This is clear from observatio ns of sys- 
tems such as M84 (jFinoguenov fc Jones! 1200 where a 
single radio outburst has disrupted the entire hot inter- 
stellar medium of the galaxy. Thus, while there is room 
to doubt the significance of heating by AGN outbursts in 
cluster cooling flows, there is a strong case that they can 
regulate the deposition of cooled gas in elliptical galaxies 
(but see IBrighenti fc Mathewj|2003|) . It would be sur- 
prising if the mechanism that regulates cooling in the 
larger systems was fundamentally different. 

If AGN feedback is the process that regulates the de- 
position of cooled gas in cooling flows, then it must be 
powered by the accretion of cooled and cooling gas onto 
the AGN. The remainder of my talk is a simpleminded 
consideration of what we should expect in that case. 

2. Bondi Acretion from Cooling Hot Gas 

It is well known that angular momentum plays a criti- 
cal role in accretion onto a black hole, i.e. an AGN. Since 
the ratio of the Bondi radius, where the influence of the 
black hole on the surrounding gas starts to be signifi- 
cant, to the Schwarzchild radius is roughly c^/s^ ^ 10^ 
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- 10^, angular momentum that is negligible at the Bondi 
radius can have a dominant effect on the subsequent 
flow. Here it is assumed that accreted gas is fed into 
a reasonably compact accretion disc, from where it is 
fed into the AGN. At the least, this means that there 
need not be a very close link between the Bondi accretion 
rate and the instantaneous accretion rate onto the AGN. 
Nevertheless, gas that is fed into the disc is assumed 
to be available to fuel the AGN. It s hould also be noted 
that, bas ed on numerical simulations. lProga fc BegelmanI 
()2003albD have cast considerable doubt on the Bondi ac- 
cretion formula. Despite this, it is employed here for 
want of a better replacemen t. It is unclear from the 
results of iProea fc BegelmanI l)2003b|) under what con- 
ditions, if any, the Bondi accretion rate may apply. In 
particular, those simulations are adiabatic and when ra- 
diative cooling becomes appreciable at the Bondi radius 
{(3^1 below) the Bondi result may be more accurate, 
and the estimate of the accreted mass would still apply. 

The ratio of specific heats in the intracluster gas is 
7 = 5/3, so that the Bondi accretion rate depends on the 
gas properties through its entropy alone. The accretion 
rate from uniform gas with density po and temperature 
To may be expressed as 



Mr = TT 
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where the black hole mass is Afh ~ lO^Afg Mq, the en- 
tropy is represented by Eq — s^pg^^^, with the sound 
speed So — 519T^^^ km s~^, the electron density corre- 
sponding to Po is n-c cmT'^ and T is kTo in keV. It should 
also be noted that for 7 = 5/3 the sonic point in the 
accretion flow occurs at R = 0, and the accretion rate is 
determined by the entropy of the gas at this point. As 
a result, the galactic environment has little effect on the 
Bondi accretion rate (and we can insert the density and 
temperature at any point in the flow into equation to 
get the accretion rate). 

Radiative cooling directly affects the entropy of the 
gas, hence the Bondi accretion rate. Expressed in terms 
of the speciflc entropy 5, the energy equation for the gas 
may be written 



dt 
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where d/dt is the convective or lagrangian time deriva- 
tive. This equation can be rewritten in terms of the 
entropy proxy, E above, as 

tit f-cool 

where the radiative cooling time is deflned as usual, 

3pkT 
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This makes the link between cooling and Bondi accretion 
very clear. Cooling reduces the entropy, as in equation 
l|3Jl, and the entropy determines the Bondi accretion rate, 
as in equation QJ . Any heat input to the gas that results 
from AGN activity would be added to the right hand side 



of equation (0), allowing us to close the feedback loop. Of 
course, the details of this process may be complex and, in 
particular, cycles of heating, accretion and cooling may 
be intermittent rather than steady. 

The behavior of cooling gas is relatively simple. Ra- 
diative losses reduce the entropy of the gas, which is 
therefore compressed by surrounding gas in order to 
maintain hydrostatic equilibrium. This causes inflow. 
The most counterintuitive feature of this flow is that, 
as a result of the compression, the temperature of the 
cooling gas is maintained at about the "virial tempera- 
ture" of the gravitational potential in which it resides. 
That is, the gas temperature at radius R is given by 
kT{R)/(pmn) ^ GM{R)/R (to within a factor of or- 
der unity), where M{R) is the gravitating mass within 
R (this only fails if the virial temperature decreases too 
rapidly with decreasing radius). 

The primary condition required to maintain the tem- 
perature of the cooling gas close to the virial temperature 
is that the gas must remain approximately hydrostatic. 
The condition for this is that the cooling time of the gas 
exceeds its sound crossing time, i.e. at radius R, that 
^cooi <^ tsc — R/s, where s is the speed of sound in the 
gas. If the cooling time falls below the sound crossing 
time, the gas will cool to low temperature, more-or-less 
in place {i.e. isochorically) , and its subsequent behavior 
is strongly time-dependent. 

For the Bondi accretion flow, what matters is the ratio 
of cooling time to sound crossing time at about the Bondi 
radius. Inside the Bondi radius, the ratio of cooling time 
to sound crossing time increases inward in the Bondi flow 
(oc T/A(T)), so that radiative cooling is unimportant 
throughout if it is unimportant at the Bondi radius. The 
ratio of cooling time to sound crossing time at the Bondi 
radius is 
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(where we have used po/'^-H.o = 4mH/3, as appropriate 
for gas that is 25 percent helium by mass). Thus, as the 
gas cools, the temperature of the gas at the Bondi radius 
remains nearly constant, at about the virial temperature 
as long as /3 > 1. If the cooling time becomes shorter 
than the sound crossing time at the Bondi radius, then 
the gas temperature at the Bondi radius will plummet. 
Formally, this causes a dramatic increase in the Bondi 
radius and accretion rate, but in practice the flow will 
become strongly time-dependent. Gas near to the Bondi 
radius will go into free-fall, and we should expect a sharp 
rise in the accretion rate onto the AGN. 

It is instructive to express the Bondi accretion rate in 
terms of /3. Using equation ||SJ) to determine the elec- 
tron density, the ratio of the Bondi accretion rate to the 
"Eddington" accretion rate is 



TUB 



MEdd /3A(To) 
where A(To) = IQ-^^A 
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23 erg cm^ s"""^, and the "Ed- 
dington" accretion rate is defined as usual by 

47rGA'fh"^HC 



^Edd 



(7) 



AGN INTERACTION 



3 



where L^dd is the Eddington luminosity and i] = O.Itj-i 
is the radiative efficiency of accretion onto the AGN. The 
significant feature to note is that the Bondi accretion rate 
from the hot gas approaches the Eddington accretion rate 
as P approaches 1. 

3. Outburst Energy 

In what follows, it is assumed that, in the absence of 
feedback from the AGN, the entropy of the gas arriving 
at the Bondi radius decreases with time due to cooling. 
However, if the Bondi accretion rate always exceeds the 
rate at which gas would cool to low temperatures in a 
simple cooling flow, then Bondi accretion can outcom- 
pete cooling and stop the entropy of the gas at the Bondi 
radius from decreasing. At least for the case of groups 
and clusters of galaxies, where the old (morphological) 
estimates of c ooling rate gr eatly exceed the Bondi accre- 
tion rate fe.g. lFabianlll994D . this cannot occur. 

While it is unclear what would bring on an AGN out- 
burst, it seems very likely that it is related to the accre- 
tion rate. For example, it has been argued for some time 
that the radiative efficiency of AGN accretion switches 
from low to high when the accretion rate exceeds a 
small fraction of t he Eddington accretion rate, M^dd 
ijAbramowicz et aLil QQS: Nar avan fc Y i 1995). This may 
also lead to a sudden increase in the mechanical energy 
output of the AGN. Alternatively, the accretion rate may 
need to reach closer to MEdd or even exceed it somewhat 
to cause an outburst. For the purpose of the argument 
that follows, the only significant issue is that an outburst 
occurs when the accretion rate reaches some multiple of 
the Eddington accretion rate that is not a lot more than 
unity. In that case, we may assume that the gas is ap- 
proximately hydrostatic at the Bondi radius up to the 
onset of the outburst. 

We can now calculate the mass of gas accreted before 
the outburst, or, more precisely, the mass of fuel available 
to power an outburst. The mass accreted in the time 
interval ti to tf is 

/■if /-Si jyi 

Ma= / MBdt^ / MBieool^, (8) 

where equation Q has been used to change variable in 
the integral, and Sj and Ef are the entropies of the ac- 
creting gas at tj and if, respectively. This change of 
variable involves an approximation, since equation Q 
applies to a single fluid element, whereas different fluid 
elements are continually being accreted. The approxi- 
mation relies on the gas entering the accretion flow be- 
ing nearly uniform. In general, it is fairly crude, but we 
are assisted by two things. First, it is likely that the ac- 
creted gas that fuels the outburst comes, initially, from a 
fairly narrow range of radii (in fractional terms) , so that 
it is likely to be reasonably uniform. Second, the sense 
of our error is to underestimate dt^ hence the total mass 
accreted. This is because the entropy is a non-decreasing 
function of the radius in a stably stratified atmosphere 
and the pressure is a decreasing function of the radius. 
Both effects mean that cooling is slower at larger radii, 
so that the decreasing entropy of an outer gas shell lags 
behind that for an inner one. 



Using the Bondi accretion rate and the cooling time 
we find that Metcooi only depends on the gas tem- 
perature. As argued above, this remains constant until 
the cooling time is comparable to the sound crossing time 
at the Bondi radius. Thus the integral gives 



Ma ~ MB.i^cool.iln:;:^. 



(9) 



While we have only a vague idea of the values to use 
for Ej and Sf, the result is quite insensitive to these. 
We therefore take Afa = X-^B,iicooi,i, where the latter 
two factors can be evaluated at any one time when the 
cooling time exceeds the sounding crossing time at the 
Bondi radius. Allowing a reasonable range of E in the 
logarithmic factor, and including some correction for the 
crude approximation involved in the change of variable 
from t to E, reasonable values of x are of the order of 10. 

Using the results from above, we get (again assuming 
25 percent helium by mass) 

48^x(GMhmH)^ ,,5 XM| y^^^ ^^^^ 



Ma 



2 X 10^ 
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using the notation from above for the numerical factors. 
If this mass is accreted onto the black hole and converted 
to mechanical energy with efficiency rymechj then the en- 
ergy released is 



-Ba = ?7mcchA4c 



3.7xl0^«^-;-^^^'erg 
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ri/2A_23 
1. This is our estimate of the 



where r/,„och = 0.1?7, 
energy available to power an outburst. 

4. Application to Hydra A 

One of the key assumptions made in the calculation 
above is that the feedback occurs in outbursts. This is, 
at least, consistent with the observation that not a ll cool- 
ing fl ows contain an active central radio source ijBurnj 
1199(1) . If the cooling flows that lack radio sources are 
in between outbursts, then this also tells us that the ev- 
idence of outbursts fades with time. Indeed, the radio 
lobe cavities, which are the main evidence of mechanical 
energy input, almost certainly rise buoyantly and disap- 
pear, either because the decreasing contrast of the rising 
cavity makes it difficult to detect, or because the bubbles 
are non-adiabatic, so that their energy is dissipated and 
the bubbles collapse. As a result, estimates of outburst 
energy based on observations of bubbles can only give 
lower limits to actual outburst energies. 

Hydra A is noteworthy for being one of the most pow- 
erful Fanaroff-Riley class I radio sources and for hav- 
ing substa ntial bu bbles associated with its radio lobes 
jMcNama ra et al.l 2000). Using the gas properties from 
iDavid et al. I (120011) ■ the central temperature is kT ~ 3 
keV and the central abundanc e \s Z ~ 1, giv i ng the 
cooling function A_23 — 1-8. ISambruiia et all Il2n0(1t) 
estimated the mass of the black hole in Hydra A as 
Mh ~ 4 X 10^ Mq. Using these numbers in equation 
l(TT)) gives an outburst energy of i?a — 2 x IQ^'^ xVra.-i 
erg. 

For the cavity created by the southwest radio lobe in 
Hydra A, the work required to inflate the cavity is pV ~ 
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2 X 10^^ erg. This must be added to the thermal energy 
in the lobe, pV/{'^ — 1), where 7 is the ratio of specific 
heats for the plasma inside the cavity, to obtain the total 
energy required to create the lobe. If the lobe plasma is 
relativistic, 7 = 4/3 and the energy required to create the 
lobe is ^pV. This quantity needs to be doubled again to 
allow for the northeast cavity, giving a total energy for 
the bubbles of about 1.6 x 10^" erg. This agrees with the 
esimate of outburst energy for the model if xHrn,-! — 8. 
Note that, if the expanding cavities drove shocks into 
the intergalactic medium, then even more energy was 
required to create them. 

Given the uncertainties, these results are broadly con- 
sistent with the outburst model described here. Clearly, 
a lot more data are required to test this model. As 
discussed above, we should expect to find evidence 
of a range of outburst energies, extending up to the 
value given by equation A more detailed model 

for the evolution of the bubbles is required to deter- 
mine the distribution of bubble energies we should ex- 
pect to observe, but we can reasonably expect bubbles 
to evolve on timescale s comparable to their rise times 
ijChurazov et alJl2002fl . Since these are comparable to 
the estimated intervals between radio outbursts, a sig- 
nificant fraction of systems are expected to show bubble 
energies comparable to the outburst energy. Of course, 
outbursts of significantly greater energy than given by 
equation could not be accounted for by this model. 

The outburst energy in Hydra A is close to the maxi- 
mum that can be accounted for using . If the accre- 
tion rate remained an order of magnitude or more smaller 
than the Bondi accretion rate throughout the time be- 
tween outbursts (Proga & Bcgelman 2003b), then the 
AGN could not accrete enough hot gas to power forma- 
tion of the bubbles in Hydra A. In that case, the AGN 
would need to be powered chiefly by accretion of cold 
gas. This result is not sensitive to the model assump- 
tions used here. 
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5. Conclusions 

Heating by AGN feedback may not be the only mech- 
anism that prevents the deposition of cool gas in cool- 
ing flows, but it almost certainly plays a significant role. 
AGN heating may be augmented by other processes, par- 
ticularly thermal conduction in clusters, but it could 
be the only heating mechanism required in smaller sys- 
tems, expecially for isolated elliptical galaxies. This 
is not to dismiss the major gaps in our understanding 
of the heatin g pro cess (e.g. Brighenti & Mathews 2003t 
iFabian et aLll200lj) . 

Because radiative cooling decreases gas entropy, di- 
rectly increasing the Bondi accretion rate, Bondi accre- 
tion of cooling gas onto a central black hole makes a 
good candidate for part of the AGN feedback process. 
Here, Bondi accretion is assumed to provide the fuel for 
AGN outbursts. The Bondi accretion rate climbs close to 
the Eddington accretion rate as the cooling time at the 
Bondi radius falls towards the sound crossing time. An 
outburst can be triggered at about this stage, perhaps 
when the radiative efficiency of accretion switches from 
low to high, or later when the luminosity reaches close to 
the Eddington limit. Accreting gas remains roughly hy- 
drostatic near to the Bondi radius as long as the cooling 
there is shorter than the sound crossing time. In that 
case the temperature of gas at the Bondi radius stays 
close to the virial temperature as the AGN accretes fuel 
for the outburst. The mass of fuel accreted is then insen- 
sitive to the details of when the outburst occurs. This 
gives an estimate for outburst energy, cauation llll which 
depends weakly on uncertain details. This outburst en- 
ergy is determined mostly by the mass of the black hole 
and the virial temperature of its host galaxy. 
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